Large-Scale Tests of the DGP Model 
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The self-accelerating braneworld model (DGP) can be tested from measurements of the expansion 
history of the universe and the formation of structure. Current constraints on the expansion history 
from supernova luminosity distances, the CMB, and the Hubble constant exclude the simplest 
flat DGP model at about 3cr. Even including spatial curvature, best-fit open DGP model is a 
marginally poorer fit to the data than flat ACDM. Moreover, its substantially different expansion 
history raises serious challenges for the model from structure formation. A dark-energy model with 
the same expansion history would predict a highly significant discrepancy with the baryon oscillation 
measurement due the high Hubble constant required. For the DGP model to satisfy this constraint 
new non-linear phenomena must correct this discrepancy. Likewise the large enhancement of CMB 
anisotropies at the lowest multipoles due to the ISW effect would require either a cut off in the initial 
power or new phenomena at the cross over scale. A prediction that is robust to both possibilities is 
that high-redshift galaxies should be substantially correlated with the CMB through the ISW effect. 
This correlation should provide a sharp test of the DGP model in the future. 
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I. INTRODUCTION 

Cosmic acceleration may arise either from a new 
form of energy density, dubbed dark energy, or from a 
modification of gravity on cosmological scales. Dvali, 
Gabadadze and Porrati (DGP) 1] have proposed a 
braneworld modification of gravity to explain accelera- 
tion (see Q for a recent review). In this model our uni- 
verse is a (3+l)-dimensional brane embedded in an infi- 
nite Minkowski bulk. The cosmological solution of this 
theory exhibits self-acceleration on the brane classically 
II- 

Various cosmological probes which can distinguish 
DGP from dark energy models (DE) have been inves- 
tigated in the literature 0, i, 0, S, @, El El GJ, El El, 
ElEBj- These come in two classes: those that test DGP 
modifications of the expansion history and those that test 
changes in how cosmological structures are formed. The 
advantage of the former class is that the predictions of 
the DGP model are well understood theoretically. The 
main disadvantage is that they currently lack the power 
to distinguish the DGP model from its dark energy coun- 
terpart the ACDM model at high significance. Moreover, 
expansion rate tests cannot distinguish between the DGP 
model and an arbitrary form of smooth dark energy that 
is constructed to mimic it in the background. Structure 
formation tests on the other hand can in principle solve 
both these problems. However due to complexities in the 
theory they currently can only be rigorously implemented 
on a rather small range of scales. 

In this paper, we examine the combination of these 
two types of tests. We employ three tests of the expan- 



sion history: the angular diameter distance to recom- 
bination from the three- year WMAP data [17j, luminos- 
ity distance to high-redshift supernovae (SNGold[l8( and 
SNLS[li|) and the local expansion rate from Hubble con- 
stant measurements. These tests do in fact strongly dis- 
favor the simplest version of the DGP model: a modi- 
fication of gravity in a spatially flat cosmology with no 
bulk cosmological constant. The best fit DGP model has 
a significant negative curvature and a high Hubble con- 
stant. Even with the addition of spatial curvature, it is 
a slightly poorer fit to the data than flat ACDM. 

The best fit DGP model yields an expansion history 
that is strikingly different from flat ACDM. The dark 
energy model that mimics this expansion history has 
negative curvature and a dark energy equation of state 
that increases sharply with redshift. In the dark energy 
context such an expansion history would be strongly ex- 
cluded by structure formation tests in weakly to strongly 
non-linear regime. These include the baryon acoustic os- 
cillations (BAO) in galaxy clustering [lfj, [20(, the abun- 
dance of galaxy clusters, and fluctuations in the Lyman- 
a forest. While these probes likely also test the DGP 
model, none of these phenomena are sufficiently well un- 
derstood in the model to provide a definitive test. Indeed 
the presence of ghost degrees of freedom in the DGP 
model around a de-Sitter background [2l|, [H [H [24[ re- 
quire that the nature of gravitational interactions change 
dramatically in the strongly-coupled non-linear regime 
for the model to remain a viable explanation of cosmic 
acceleration 
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Under this assumption, structure formation in DGP is 
currently well understood on scales between the Hubble 
scale and the scale radius of a typical dark matter halo. 
Large scale CMB anisotropy is substantially enhanced by 
the Integrated Sachs- Wolfe (ISW) effect. At the lowest 
multipoles, this excess poses a significant challenge for 
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the model. Possible resolutions of this problem still leave 
a robust prediction and sharp test of the DGP model: 
high redshift galaxies should be substantially correlated 
with the CMB. 

II. EXPANSION HISTORY 

A. Supernovae, CMB and Ho 



We implement the SN constraints from the "gold" SN 
data set of and the SNLS set of [H]. For the CMB, 
we fix the distance to recombination at z\ ss = 1088^ 
through the measurement of the acoustic peak scale 
I A = 302^4 an d its length calibration through the mat- 
ter density fl m h 2 = 0. 1268±g;gg|§ 1293. Note that the SN 
are a relative distance indicator and hence constrain D 
whereas the CMB is an absolute distance measure which 
constrains D/h. 



The expansion history of a DGP model is deter- 
mined by the usual Friedman-Robertson- Walker param- 
eters and the crossover distance. The crossover distance 
is defined as the ratio of 5-dimensional to 4-dimensional 
Planck mass scales 
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It is a free parameter of the theory. If r c is close to the 
current horizon scale, the acceleration of cosmic expan- 
sion is replicated without dark energy. Gravity on scales 
much less than r c does not penetrate into the bulk signifi- 
cantly, and exhibits the usual Newtonian potential falling 
as the inverse distance. In the DGP model, as the Hubble 
scale approaches r c , the universe enters an accelerating, 
asymptotically de-Sitter phase despite the absence of a 
cosmological constant. 

The main evidence for cosmic acceleration is geometri- 
cal and based on the luminosity distance to high-redshift 
supernovae (SN) [l8|, [2(| [2?], [28[ and the angular diame- 
ter distance to the recombination epoch of the CMB (e.g. 
(29|). Both distances are given as usual by the (comov- 
ing) angular diameter distance in units of H^ 1 
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FIG. 1: The A\ 2 between the best fit flat and open DGP versus 
that of flat ACDM model. SNGold supernova (SN) data set is used 
in the top panel and SNLS SN data set is used in the bottom panel. 
The DGP model requires curvature and a high Hubble constant. 
With the addition of Key Project (KP) direct Hubble constant 
measurements open DGP is a marginally poorer fit to the data 
than flat ACDM. 



where r(z) is 
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E(z) is the Hubble parameter in units of the Hubble con- 
stant. The DGP model modifies the Friedman equation 
or the relationship between E(z) and the energy densities 
of matter and radiation in units of the critical density 
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where Hq = 2997.9ft. Mpc 1 . Here the scale factor a 
1/(1 + z), a c (=l/4r c 2 i? 2 ) is given by 



4(1 -fit) 
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Like flat ACDM (fACDM), flat DGP has a single pa- 
rameter Q m that is over-constrained by the distances to 
recombination, fl m h 2 , and high-redshift supernovae. Un- 
like fACDM, flat DGP is not a good fit to these data. 
In DGP cosmologies, the SN favor a lower VL m than in 
fACDM resulting in a shorter angular diameter distance 
to the CMB than is measured by WMAP. For the SNGold 
data set the excess x 2 over fACDM is A\ 2 ss 4 while for 
SNLS the excess is A% 2 « 9 (see Fig.[TJ(. 

A negatively curved geometry elongates the CMB dis- 
tance and allows a low f2 m in agreement with the SN dis- 
tances [l6l. [30|. Even so, the CMB constraint on fl m h 2 
favors a high Hubble constant of h = 0.8 (SNGold) and 
h 0.9 (SNLS) compared with observations from the 
HST Key Project (KP) h = 0.72 ± 0.08 fACDM 
on the other hand, despite being highly over-constrained 
satisfies the observations with h = 0.73. 

With the addition of the KP Hq constraint, the best 
fit open DGP (oDGP) model has h = 0.76 (SNGold) and 
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h = 0.82 (SNLS) but is a marginally poorer fit to the 
data than fACDM. The excess over fACDM of the latter 
is A% 2 ~ 2.5 with one fewer degree of freedom due to 
the fit to spatial curvature. The other parameters of this 
best fit model are Q m = 0.18, Qk = 0.036. 

If the Hubble constant measurements can be improved 
by a factor of 2 in the future, then the oDGP and fACDM 
will be clearly distinguished from each other. They are 
currently both acceptable fits to these data. We will 
hereafter take the best fit model to the WMAP, SNLS 
and KP data sets as the oDGP model for further tests. 



B. DGP vs. Dark Energy 

The main expansion history differences between the 
best fit oDGP model and fACDM can be exposed by 
finding the scalar field or QCDM model that matches 
the oDGP distance measures exactly. With the same 
curvature, il m h 2 and Hq, that QCDM model has a time- 
dependent equation of state 
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This equation of state starts at w = —1/3 during radia- 
tion domination, becomes w = —0.5 during matter domi- 
nation and decays to w — —0.85 at the present time. The 
value of w(a) is substantially higher than the w = — 1 of 
fACDM. The implied increase in the dark energy density 
causes both the larger Hubble constant and a decrement 
in D(z) at high redshift. This decrement is compensated 
by negative curvature to match the CMB distance. 

The expansion history of this QCDM model is strik- 
ingly different from fACDM despite the fact that they 
both satisfy CMB and SN constraints. Both the in- 
creased dark energy density and the negative spatial cur- 
vature contribute to changes in structure formation that 
would violate even current observations. One therefore 
expects that structure formation tests should provide a 
strong constraint on the oDGP model. It is important to 
emphasize that it is not necessary to distinguish between 
oDGP and QCDM based on difference in the growth of 
structure. The comparison is between oDGP and fACDM 
because the QCDM model is not viable. We therefore use 
the degenerate QCDM model as a guide to understand 
the qualitative predictions of oDGP. The main task of 
testing the expansion history consequences for oDGP is 
to find a structure formation probe that can be robustly 
calculated in the DGP model. 



III. STRUCTURE FORMATION 

A. Baryon acoustic oscillations 

As discussed in the previous section, the largest ex- 
pansion history discrepancy between oDGP and fACDM 



lies in the current expansion rate or Hubble constant. In 
the distance-degenerate QCDM model, the large Hub- 
ble constant would be in strong conflict with the mea- 
sured baryon acoustic oscillations in the SDSS galaxy 
correlations. Baryon acoustic oscillations appear near 
the scale where density perturbations become non-linear 
and hence cannot be considered a purely geometric test. 
Detailed predictions in the oDGP model below this scale 
require the solution of the A^-body problem with a non- 
linear Poisson equation. It is nonetheless useful to quan- 
tify the extent of the problem in the QCDM model since 
it determines the level at which nonlinear oDGP cor- 
rections would have to alter linear theory predictions to 
satisfy the data (l6| . 

The baryon oscillations are imprinted at a fixed phys- 
ical scale by the CMB at recombination. Under ordi- 
nary gravity, they yield a nearly geometric test since non- 
linearities are insufficient to change the physical scale of 
features substantially. The SDSS luminous red galaxy 
(LRG) survey measures the features between 0.16 < 
z < 0.47. At these redshifts, this distance is sensitive 
mainly to the Hubble constant [32]. More specifically, 
it can be converted to a relative distance measure be- 
tween the effective LRG redshift and recombination of 
D(z = 0.35)/L»(z lss ) = 0.0979 ± 0.0036. Recall that the 
CMB measures D(z\ ss )/ h. Hence in the high Hubble con- 
stant QCDM model this number is too low, specifically 
0.0818 leading to a 4.5er discrepancy. This discrepancy 
can be slightly ameliorated by refitting the models to all 
data including the BAO data 

If oDGP shared the same non-linear effects with 
QCDM, the BAO measurement would exclude it as well. 
In the DGP model, there is a key distance scale, r* , above 
which the linearized perturbations experience the weak- 
brane phase @, [H, Hj, [H, [36[ . In the weak-brane regime 
r c r ^> r», gravity is modified as a scalar-tensor the- 
ory. Below , gravity is described by the standard Ein- 
stein theory. For a point mass M, the transition scale 



r* = (r 2 c r g )^ 3 , (r g - 2G N M). 



(7) 



Since even the linear density field is constructed from the 
spatial average of collapsed dark matter halos under any 
hierarchical structure formation model, the critical scale 
is associated with the of the dark matter halos. 

For definiteness let us take the other cosmological pa- 
rameters as in the best fit QCDM model: VL m h? — 0.122, 
ilbh 2 = 0.0224, r = 0.11 with a power law spectrum of 
adiabatic fluctuations with n s = 0.958, 6{ = 4.5 x 10 -5 
at fc„orm = 0.05 Mpc -1 . In this QCDM model, the 
LRG galaxy number densities correspond to halos of 
M* « lO 13 /i _1 Af corresponding to r* s» 3 Mpc. This 
scale is only marginally smaller than the baryon oscil- 
lation scale. Furthermore, without a determination of 
how dark matter halos form under DGP and whether 
the observed galaxies indeed reside in 10 13 /i _1 M Q halos 
it remains unclear whether weak-brane phase predictions 
can be applied to the BAO measurements. 
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B. Linear Growth Rate 

The same expansion history discrepancy between 
QCDM and fACDM which raises h also causes both an 
earlier decay of the growth of structure and a larger to- 
tal change in the growth to z = Q. For example, in the 
QCDM model <t 8 = 0.6 and fi ro = 0.18. These values 
are in conflict with the abundance of local clusters. Un- 
fortunately, this QCDM prediction again relies on the 
non-linear aspects of structure formation and hence can- 
not yet be applied to oDGP. In order to test oDGP we 
must find a probe on scales r c r ^> r* . 

In this regime, the linear theory predictions of DGP 
are now well understood. In general, the equations of 
motion of linear perturbations in DGP are not closed on 
the brane. They require knowledge of the gradient of 
metric perturbations into the bulk as quantified by the 
so-called master variable [37j ■ This in turn depends on 
boundary conditions in the bulk. Koyama and Maartcns 
[HI solved this system in the small scale limit k / aH 3> 1 
where the expansion rate can be ignored compared to 
the master variable frequency u> ~ k. In this quasi-static 
(QS) limit, the Newtonian metric perturbations become 
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during the transition from matter domination to acceler- 
ated expansion but the iteration yields stable solutions. 
We call this the dynamical scaling (DS) solution. 

Under these assumptions the quasistatic solution is 
stably approached by k = lOaH and hence is valid to 
better than 2% at k > 0.01 Mpc -1 for 2 = 0. Com- 
bined with the condition from r*, a robust test of the 
quasi-static predictions of oDGP would involve fluctua- 
tions between 0.01 < ^(Mpc" 1 ) < 0.1. 



where <&_ = ($- *)/2 and $+ = ($ + *)/2 and 
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Here A,„ is the comoving matter density perturbation 
and for k/aH ^> 1 satisfies 



k 2 

A m + 2HA m = -(^ 



(10) 



These relations agree with the earlier conclusions of Lue 
obtained in a different manner. 

This QS solution implies a 5% to 10% extra decay of 
$_ in the oDGP compared with the QCDM model with 
the same expansion history. We show the evolution of 
$_ for each model in Fig. [5] The main effect however is 
not a difference between oDGP and QCDM but between 
either and fACDM. In either case, <E>_ begins to decay at 
high-z such that there is as much decay at z > 1 as there 
is in fACDM at z < 1. 

In order to find an appropriate test of this prediction, 
we need to determine the validity of the quasistatic so- 
lution on large scales approaching r c . In a companion 
paper {38j . we solved the dynamical equations by mak- 
ing a starting assumption that outside the horizon the 
evolution is scale free and that perturbations in the bulk 
vanish at the causal horizon. We then iterate the solution 
through the dynamical equations until they reach conver- 
gence. The starting scaling assumption is mildly violated 



FIG. 2: The growth factor of 3>_ of fACDM and oDGP under 
the quasi-static (QS) assumption normalized to <3>_(z 3> 1) = 1. 
The model parameters are chosen as the best fit to the CMB (3yr- 
WMAP), SN (SNLS), and H (KP). QCDM is the scalar field dark 
energy model which has the same expansion history as oDGP. The 
growth rate is plotted on subhorizon scales where the dark energy 
is smooth and oDGP follows the QS assumption. 



C. ISW Effect 

The decay in the gravitational potential <&_ at k > 
0.005 Mpc -1 causes an enhancement of the large-angle 
anisotropy in the CMB through the ISW effect. 

The angular power spectrum of the ISW effect is given 

by 



CP = 4*|^[I?(fc)] : 



,fc 3 P t 
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where P$„$_(fc,0) is the power spectrum of $_ at the 
present time and the kernel if is 



(12) 
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The window function is given by 

2 <9$_ 



W\k,z) = 



$_ dz 



(13) 



We have have here assumed that the spatial curvature is 
small enough that the spherical Bessel function ji accu- 
rately represents the radial harmonics. 
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FIG. 4: Top panel: the spatial wavenumbers corresponding to the 
ISW effect under the QS approximation for multipoles I = 2, 20. 
For I = 20 the scales associated with the effect are l/r c <Cfc'Cl/»"* 
and even the quadrupole receives much of its contribution from 
l/r c -C k. Bottom panel: the fractional difference between the 
dynamical scaling (DS) and quasi-static (QS) solutions for DGP at 
the relevant scales for I = 2, 20. For I = 20, the QS approximation 
is accurate to the several percent level. 



FIG. 3: The CMB temperature anisotropy power spectra of the 
various models. Top panel: oDGP under the quasi-static (QS) ap- 
proximation versus QCDM and fACDM for the total power (upper 
curves) and the ISW effect (lower curves). Bottom panel: oDGP 
(QS) compared with the dynamical scaling (DS) solution. The DS 
solution predicts a further enhancement of the ISW effect such that 
I < 10 strongly violates the WMAP measurements (points). The 
errors are mainly associated with cosmic variance at these multi- 
poles and the 68% and 95% CL bands per I have been attached to 
the DS model. 

Fig. [3] shows the difference between predictions for Cj 1 
from oDGP under the quasistatic approximation com- 
pared with fACDM and QCDM. We also show the to- 
tal temperature power spectrum C; TT . The enhance- 
ment of the low-multipole power spectrum directly re- 
flects the high-redshift decay of the gravitational poten- 
tial in oDGP and QCDM. In QCDM the effect at the low- 
est multipoles is smaller since a scalar field has a Jeans 
scale at the horizon. Beyond this scale, dark energy per- 
turbations slow the decay of the potential. 

There are analogous corrections for oDGP on large 
scales but they have the opposite sign. Under the DS ap- 
proximation the decay rate of the actually increases 
at the largest scales (see Fig. H]) leading to even larger 
fluctuations. In [38| . we show that this approximation 
should at least give the right sign of the effect. The QS 
solution should thus be taken as a lower bound on the 
ISW effect. Hence the discrepancy between the oDGP 



QS predictions and the measurements at I = 2 — 5 in the 
3year WMAP data are a challenge to the model that will 
require either a cut off in the initial spectrum or other 
new phenomena near r c to overcome. 

At I ~ 20 the ISW effect comes from physical scales 
that are within the limit of the quasistatic approxima- 
tion as well as well below r c . The ISW effect here is a 
small fraction of the total anisotropy and is best isolated 
through cross correlation with tracers of the gravitational 
potential, e.g. galaxies. Moreover the cross correlation 
can isolate the redshift history of the potential decay and 
expose a key prediction of the oDGP structure formation: 
an early decay of the potential. 



IV. GALAXY-ISW CORRELATION TEST 

The oDGP model predicts a stronger and earlier decay 
of the gravitational potential than the fACDM model. 
By cross-correlating galaxies of different redshifts with 
the CMB, one can in principle reconstruct the redshift 
history of the potential decay. Furthermore, the cross- 
correlation arises from the well understood quasi-static 
(QS) regime of oDGP. 

The cross-power spectrum of the CMB and a set of 
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FIG. 5: The galaxy-ISW cross-correlation coefficient R\ in each galaxy bin from z — to z = 3. Solid curves denotes fACDM 
and dash curves denotes oDGP. Note the much larger correlation at high z in oDGP. 



galaxies gi is given by a generalization of Eq. (jlip 
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where the galaxy kernel I 9i is 



1/2 



.(15) 



Under the QS approximation, the window function be- 
comes 



W 9i (k,z) 



2 k 2 m{z)bi{z) 
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where rii(z) is the redshift distribution of the galaxies 
normalized to J dzrii — 1 and bi(z) is the galaxy bias. 

For definiteness, we assume that the galaxy sets come 
from a net galaxy distributions of 



n g (z) oc ^e-^ 1 ' 5 ) 2 , 



(17) 



where the normalization is given by the LSST expecta- 
tion [42| of 35 galaxies per arcmin 2 . For the subsets of 
galaxies, we assume that this total distribution is sep- 
arated by photometric redshifts which have a Gaussian 
error distribution with rms a(z) — 0.03(1 + z) (see [39[ 
for details). The redshift distributions are then given by 



A- 



crfc 



-1-2 



(t 



— erfc 



/ Zj - z 

\V2a(z) 



where erfc is the complementary error function and Ai is 
determined by the normalization constraint. 

Our galaxy bias bi(z) is determined by a halo model 
and ranges from 1.4 at z ~ to 3.8 at z ~ 3. How- 
ever with our narrow binning, the bias is nearly constant 



across the bin and so its value can be empirically de- 
termined through comparing the auto and cross power 
spectra once the data are in hand. In particular, the 
cross-correlation coefficient R\ in harmonic space 



R) 



Cf* 



(18) 



is independent of the bias and is a robust measure of the 
galaxy-ISW correlation. 

In Fig. El we show R\ for oDGP and fACDM. The main 
difference appears in the high-redshift bins and reflects 
the early decay of the potential in oDGP. The correla- 
tion of high-redshift galaxies with the CMB is therefore 
a sharp test of the oDGP model. 

To determine whether these correlations are poten- 
tially observable, we estimate the sampling and noise er- 
rors in the cross-correlation measurement. First let us 
define the total galaxy power spectra Cf' 9 '' as 



Nt 9i 



(19) 



Here the shot noise Nf i9i = 1 jn % A where n l A is the angular 
number density of galaxies in the bin. With our number 
densities the variance is nearly sample limited in all of 
the redshift bins. 

Likewise the total power in the CMB measurements is 
given by 



TT 



c, 



TT 



(20) 



where the instrumental noise -/Vj TT is A^e'^ 1 ^ 2 / 81112 . 
With either WMAP or Planck satellite assumptions for 
the noise At and the beam a the temperature map is 
sample variance limited to good approximation. 

The error on cross-correlation can be calculated under 
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the Gaussian assumption as 



ACfil -V(2l+l)f sky 



-,1/2 



,(21) 



where / s k y is the sky fraction covered by the galaxy sur- 
vey. 

Statistical Properties of Galaxy Correlations 







z 


0.2 1 0.6 1 1.0 1 1.4 1 1.8 1 2.2 1 2.6 1 3.0 




2.5 1 5.0 1 5.9 1 6.0 1 5.7 1 5.3 1 4.9 1 4.5 




12 | 18 | 23 | 26 | 29 | 31 | 33 | 35 


noDGP / nfACDM 

H ho% / H ho% 


0.83 1.1 1.6 2.2 2.8 3.6 4.4 5.3 







TABLE I: Rows: mean redshift, total signal-to-nose of corre- 
lation, multipole I below which 50% of the (S l /N 1 ) 2 arises, 
ratio of correlation between oDGP and fACDM at that mul- 
tipole. 



The signal to noise in each bin is given by 




(22) 



and these values are given in Tab. HI In the sample vari- 
ance limit approached by our fiducial assumptions 



m 2 



(23) 



Even at high redshift z ~ 2 — 3 the signal to noise allows 
a 15% — 20% measurement of the correlation. Such a 
measurement would be more than sufficient to distinguish 
the factor of a few difference between the oDGP and 
fACDM models. 

Most of the detectable signal comes from the regime 
where the QS approximation is appropriate. At low mul- 
tipoles I < 10, the cosmic variance of both the galaxy 
and temperature fields dominate [40| . By I = 100 the 
sample variance of the acoustic peaks dominates. From 
the signal to noise sum in Eq. (|22[) . one can define Z 50 % as 
the multipole at which the (S l /N 1 ) 2 reaches half its total 
value. This scale is also given in Tab. U and lies between 
I = 20 — 40 for most of the bins. At this scale the ratio of 
correlation coefficients in oDGP and fACDM approaches 
a factor of 6 at z = 3. High redshift galaxies can thus test 
the quasi-static oDGP predictions with high significance. 



supernovae and the Hubble constant. Unlike a model 
with a cosmological constant as the dark energy, the DGP 
model requires spatial curvature to satisfy the CMB and 
supernovae constraint and even so mildly exceeds bounds 
on the Hubble constant. 

In a dark energy model, the DGP expansion history 
would be ruled out by structure formation tests. These 
are serious but perhaps not insurmountable challenges 
for the model. Among the most powerful are the baryon 
oscillations and the large angle CMB anisotropy from the 
ISW effect. In the DGP model, these two phenomena lie 
close to the two critical scales: r* where gravity tran- 
sitions from obeying the usual Einstein equations to a 
scalar-tensor or weak brane phase and r c when gravity 
exits the weak brane phase and becomes fully 5 dimen- 
sional. At the transition, the theoretical framework 
to calculate structure formation is currently incomplete. 
At the r c scale, changes in the initial power spectrum or 
new gravitational physics may ameliorate the problem. 

The enhancement of the ISW effect from the decay 
of the gravitational potential extends to intermediate 
scales where the theoretical calculations are robust. The 
main difference between the DGP and the flat cosmolog- 
ical constant expansion history is the prediction that the 
gravitational begins its decay at high redshift. A sharp 
test of the DGP scenario would be to isolate this de- 
cay at high redshift by cross correlating the CMB with 
a high-rcdshift galaxy population. The cross correlation 
coefficient is also relatively robust to changes in the initial 
amplitude of fluctuations. We show that a galaxy survey 
that covers a substantial fraction of the sky at z > 1 can 
detect this correlation at high significance. Such surveys 
are currently being planned for the measurement of high- 
redshift baryon oscillations. Alternately, radio sources, 
X-ray sources and quasars provide other high-redshift 
populations for cross correlation. 

A possible source of confusion for this test is the ap- 
parent number density fluctuation induced by magnifi- 
cation from gravitational lensing by the same structures 
associated with the ISW effect [41|. Magnification both 
lowers the number density for a fixed population behind 
a lens and also raises it for a fixed flux limit by bring- 
ing faint galaxies into the sample. This effect may be 
separated from a true spatial correlation by examining 
different populations of objects where the competition 
between these two effects differ and the cross-correlation 
between high and low redshift galaxy populations. A full 
examination of this issue is beyond the scope of this paper 
but it is likely that the galaxy-ISW cross correlation will 
become a sharp test of the DGP braneworld acceleration 
model in the future. 



V. DISCUSSION 

The self-accelerating DGP braneworld model is chal- 
lenged by expansion history constraints from the CMB, 
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